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Abstract 

We reconsider the problem of the static thermal roughening of an elastic manifold 
at the critical dimension d = 2 in a periodic potential, using a perturbative Functional 
Renormalization Group approach. Our aim is to describe the effective potential seen 
by the manifold below the roughening temperature on large length scales. We obtain 
analytically a flow equation for the potential and surface tension of the manifold, valid 
at all temperatures. On a length scale L, the renormalized potential is made up of a 
succession of quasi parabolic wells, matching onto one another in a singular region of 
width ~ L~^l^ for large L. We also obtain numerically the step energy as a function of 
temperature, and relate our results to the existing experimental data on ^He. Finally, 
we sketch the scenario expected for an arbitrary dimension d <2 and examine the case 
of a non local elasticity which is realized physically for the contact line. 
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1 Introduction 

The roughening transition has been studied in great detail, both theoretically and experi- 
mentally 0]. Direct analogies with the (two dimensional) XF- model or the Coulomb gas 
furthermore make this problem particularly enticing 0. More recently, the role of disorder 
on the roughening transition or on the properties of the XY model, has attracted consider- 
able interest 0, ^ ^ . In particular, replica calculations and Functional Renormalization 
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Group (frg) methods have been apphed to this problem, with sometimes conflicting results 
[0. In this paper, we wish to reconsider the problem of the roughening transition in the ab- 
sence of disorder, from a frg point of view, where the flow is not a priori projected onto the 
first harmonic of the periodic potential. Within a local renormalization scheme, we establish 
exact equations for the evolution of the full periodic potential V{ip), and the surface tension 
7 with the length scale L = e^, which we analyze both numerically and analytically, in the 
low temperature phase. If we start with a sinusoidal periodic potential, the shape of the fixed 
point potential V*{ip) evolves to a nearly parabolic shape with matching points becoming 
more and more singular as the length scale increases. The nature of the singularity is inves- 
tigated in detail close to the fixed point, that is for small values of the rescaled temperature 
_ T 

T = -r where A is the periodicity of the potential and 7 the elastic stiffness. We find 

27r7A"' 

that the width A^p of the singular region scales as L"^^*^-^)/^, where g{T) governs the scaling 
of the surface stiffness with the length scale according to 'y{L) ~ L^^'^\ The exponent g(T) 
tends towards 2 with negative corrections which we calculate, when T goes to (i.e. for 
L — > 00). 

The paper is organized as follows. In section 2, we introduce the model: we outline the 
calculations involved and discuss the differences with the approach of Nozieres and Gallet 
(ng), and briefly examine the problem for d < 2. We then explain in section 3, by a 
mean field argument the origin of the singularity that develops during the renormalization 
flow. In section 4, we present a scaling form for the renormalized potential, around its 
maxima and close to the fixed point, which accounts for the nature of the singularity. Using 
our renormalization group flow, we compute in section 5 the step energy as a function of 
temperature. Finally, in section 6, we look at the case of a contact line in a periodic potential, 
as this is a physical realization of a non local elastic stiffness. 



2 Model and functional renormalization group 

We consider an elastic interface whose height fluctuations are described by a profile $(a;), 
where a: is a d-dimensional vector, in the presence of a deterministic periodic potential V. 
Supposing that the slope of the interface is everywhere small, the energy of the system is: 

nm = \jd''x {y^x)f + j (fx ^(^) (1) 

where 7 is the elastic stiffness and A the periodicity of the potential. In the absence of periodic 
potential, the height fluctuations of the surface on a length scale L scale as L'^~'^. For d > 2, 
the interface is therefore always flat. For the critical dimension d = 2, the interface is 
rough only if the temperature exceeds a certain critical temperature Tr. When the potential 
V is harmonic, this model is the continuous version of the Sine-Gordon model. NG have 
studied the statics of this problem using a two-parameter renormalization group scheme, 
and have written flow equations for 7(//) and the amplitude Vo{L) of the periodic potential. 
They suppose that during the flow, Vo remains small compared with the temperature and 
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neglect all higher harmonics of the potential. Correspondingly, within this procedure, the 
renormalization scheme ceases to be valid when Vg becomes of the order of the temperature. 
In the low temperature 'flat' phase, this occurs after a finite renormalization since Vo grows 
with distance. 

In our calculation, we consider a general periodic function with the only constraint that 
it should be sufficiently smooth (we shall explain this more quantitatively in the following). 
Since we re-sum the whole perturbation expansion in Vg/T, there is however no constraint 
on the amplitude of the potential, and the renormalization procedure can be carried on any 
length scale without interruption. The relevant coupling constant appears to be Vg/i rather 
than Vo/T. During the renormalization flow, we keep track of the whole function V{(f) instead 
of projecting onto the first harmonic, so that we have a more quantitative knowledge of the 
behaviour of the potential for low temperatures. 

Technically, we proceed by considering the partition function: 

J d^x) e-^^[*(^)l (2) 

We perform the renormalization procedure by splitting the field $ into a slowly-varying and 
a rapidly- varying part as: 

=$<(x) + $>(x) (3) 

The Fourier modes k of are such that < |A;| < |A|/s, and those of such that 
|A|/s < \k\ < |A|, where s = e'^^, |A| being a high momentum cut-off, of the order of 1/a, 
where a is the lattice spacing. We integrate over the fast modes in the partition function and 
retain only the terms that renormalize the gradient term and the potential term. The other 
terms that are generated are discarded as irrelevant. Within this renormalization scheme, 

our calculation is exact. After some algebra detailed in appendix A, we obtain a set of flow 

— V - T 
equations for d — 2, ior V = , „. , and T — — ■, the rescaled potential and temperature 

7A^|A|^ ZTTjX^ 

(note that V and T are dimensionless) : 



dT = 



where g is given by: 



« (i+r(»,))° ' {i+r{^)) 



g = 47r I d(p — + - d(p -j (5) 
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These equations call for some comments. 



The relevant perturbative parameter appears to be V, rather than V/T. In the limit 
<^ 1, and in the case where the potential is purely harmonic (i.e. V{ip) = Vo cos(27r(y9)), 
the RG equations read: 



dUo / ttT \ 



di V 7A' 







Un 

(6) 



d£ \ 7A2 J 72 A4 

where Uo = fo/|Ap. The first equation is trivial and identical to the one in NG, and im- 
mediately leads to the value of the roughening transition temperature: Tr = 27ooA^/7r, 
where 700 is the renormalized value of 7. The second is close to, but different from the 
one obtained in the particular renormalization scheme used by NG: near the critical 
temperature Tr, the coefficient between parenthesis is equal to 4 in our case and to 0.4 
according to NG. 

The renormalization of the surface tension, as measured hj g, is always positive. One 
can check that, as has been pointed out by NG, if the initial potential is parabolic (i.e. 
V{(f)) = fo0^), then the coefficient g vanishes identically, and there is no renormalization 
of the surface tension. This is indeed expected since in this (quadratic) case, all modes 
are decoupled. 

The flow equations only make sense if l^" > —1. We have checked numerically that 
if this condition is satisfied at the beginning, it prevails throughout the flow. On 
the other hand, if the initial potential is so steep that this condition is violated, the 
perturbative calculation is meaningless. This comes from the fact that metastable 
states, where the surface zig-zags between nearby minima of the potential, appear at 
the smallest length scales. In this respect, it is useful to note that the last term of the 
flow equation on V comes from the integration of the Gaussian fluctuations of the fast 
field around the slow field. The condition V > — 1 is a stability condition for these fast 
modes. If the unrenormalized potential is harmonic (i.e. V{(p) = VoCos{27np)), and the 

f o / 2-71 

7 

V 

which simplifies to — < 1 in the case where A = a. If the initial value of the potential 

7 

is too large, one actually expects the transition to become first order (but see |^). 
Actually, a variational calculation indeed predicts the transition to become first order 

when — > 1 [|10|. 



unrenormalized surface tension given by 70, then this condition reads ~ j ^ 



7 

The fundamentally new term in the above equation is the second term, proportional to 
V , and independent of temperature. This term leads to the appearance of singularities 
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in the flow equation: up to second order in V, this equation is close to the Burgers' 
equation (see below) for which it is well known that shocks develop in time. The fact 
that this term survives even in the zero temperature limit is at first sight strange, 
since one could argue that for T = 0"*", there are no longer any thermal fluctuations, 
and thus no renormalization. This argument is not correct because we are computing 
a partition function, thereby implicitly assuming that the infinite time limit is taken 
before the zero temperature limit. Such a non trivial renormalization has also been 



found in the context of pinned manifolds [11, 12|, and can be understood very simply 



using a mean-field approximation, which we detail in the next section. 

For completeness let us consider the case 1 < d < 2. By simple scaling arguments, 
we can see that for d < 1 the interface is always rough. For 1 < < 2, we have a 
roughening transition between a flat phase and a rough phase. In this case, we allow 
$ to renormalize and suppose that A renormalizes in the same way according to: 

— V - 

The other flow equations in terms of the rescaled parameters V = and T = 
(with Kd = Sd/{2nY where Sd is the (i- dimensional sphere) now read: 



7A 



f = (d-«,-2C)V'-.-^ + ^ln(l + V" 



dj 
d£ 

dT 
~dl 

where gd is given by: 



{Qd + d- 2)7 

-{g, + d-2 + 20f 



2 

9d = ^9, (9) 

with g given by equation (|^) above. These equations have a non trivial fixed point 

for gd = 2 — d and ( = 0. This corresponds to a rescaled temperature Tr and a 

renormalized rescaled potential V such that equations (P) and (^ are satisfied. For 

1 < c? < 2, we obtain Tr numerically by proceeding as follows: we self-consistently 

— dV 
solve the differential equation on V obtained by putting = for different fixed 

rescaled temperatures T, imposing that g is given by equation (|^). This enables us to 
plot g as a. function of T. The rescaled temperature corresponding to the transition 
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Figure 1: g a.s a function of the rescaled temperature T for d=3/2. The point where g = 
2 — d = 1/2 determines the transition temperature. 

temperature is such that ga = 2 — d. In the figure (|1]), we have plotted the result for 
d = 3/2, with 7o = 1. In that case, T ~ 1.2. 



3 Mean field analysis and effective potential 

In this section, we show on a simplified mean field version of the model how the non linear 
term V'"^ arises in the flow equation of the potential. Using a discrete formulation of the 
problem and replacing the local elasticity modeled by the surface tension term by a coupling 
to all neighbours, we can rewrite the energy as: 

N N 

Hmf{{<^h, = ^ E(*^ -W + aJ2 ^(*^) (10) 

i i 

where $ is the center of mass of the system, a the lattice spacing and L = Na. Implementing 
the constraint $ = 1/A^^^ $j by means of a Lagrange multiplier in the partition function, 
we have: 
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= dr] l[d<l>i e i (11) 



which can also be expressed as: 



2 

N log zi^.rf) + A^;^ ^ 



Z[$] = j dT]e 2/37 (12) 

where 



2a V (3j/ = z(^<i, + (13) 

We are left with a simpler problem since we now have a one-body problem. We introduce 
the auxiliary partition function -2_r(\E', t) defined as: 



/37 ($ - $) 



2 



, /37 f -—- --(3aV(^) 



27rar _ 

Up to a multiphcative constant, one has z[(^ + — ) = Z/jf $ + — , r = 1 ), where 2:r(\&, r) 

V P7/ V p7 / 

verifies the diffusion equation: 

dzR _ a d'^ZR 

dr ~ ^ > 

with an initial condition given by: 

^^(^,r = 0) = e-/5«^W (16) 
Defining now the effective pinning potential Vr as: 

<'^«(* + ^.-)=-riog--«(* + |^,r) (17) 
we can then easily show that Vr is the Hopf-Cole solution of the non linear Burgers' equation 
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(18) 



d{Ta) 7 ^ V 

where the temperature independent non linear term indeed appears. 

It is easy to show that when ^ 00 or T ^ 0, the original partition function can be 
solved by a saddle point method, leading after a change of variables to an effective potential 
per unit length: 
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Veff{^) = Vr{u, r = l)-^^{^-uf (19) 

where u is given by: 

f (^-^) = ^(^,^ = 1) (20) 
la oif ip=u 

Now, by changing Vr to — Vr, one can see that —Veg can also be written as the solution of a 
Burgers' equation. It is known from results on the Burgers' equation, that, with 'time' r, the 
effective potential Vr develops shocks, smoothed out at finite temperature, between which it 
has a parabolic shape. The appearance of singularities is due to the non linear term in the 
partial differential equation, which indeed survives in the limit T = 0. It is interesting to see 
how this 'toy' renormalization group captures some important features of the full scheme, 
such as the one shown above for a non disordered potential. 

4 Analysis for small T 

In this section, we go back to the model introduced in section 2 and analyze the nature of V 
close to the low temperature fixed point, that is for small values of the rescaled temperature 
T. Since 5^ > in the low temperature phase, this corresponds to the large scale structure of 
the renormalized potential for all temperatures T < Tr. 

Expanding V around one of its minima as V{ip) = Vm + ~ f*)"^: replacing V in 
the flow equation(l), we have 

(21) 

d/Ki , , Ki 

— = (2-g)K-2TT 

One can actually check that a parabolic shape for V is exactly preserved by the renormal- 
ization flow. However, since the potential has to be periodic, these parabolas should match 
periodically around each maximum that is for ^p — ^p* = 0,1, 2,.... The region of the maximum 
is therefore expected to be singular. To investigate the nature of the renormalized periodic 
potential around its maximum value, we will thus make a scaling ansatz on V for small T. 
For our perturbative calculation to be valid, we expect V (0) to be > —1. Now since we 
expect a singularity to develop as T goes to 0, it is probable (and actually self-consistently 
checked) that V (0) should tend towards —1. As T goes to 0, we thus make the scaling 
ansatz: 

l+F"(^)=TVm (22) 



where JF'(O) > 0. This means that the width of the singular region behaves as Aip ~ T". 
Hence, in the scaling region: 
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F'M = -^ + T^^"^f|,) (23) 



with JF(0) = to ensure that = is a maximum of V. Integrating once more the above 
equation, one finds: 

F(^)=FM-^+T^«"a(|,) (24) 

with Q' = T . Replacing this last equation in the fiow equation for V, we obtain: 

^ = (2-,)FM + |logT^ (25) 

Suppose that equations ( PT| ) have a fixed point as T goes to zero, and that close to the 
fixed point one can neglect the left hand side of these equations. This leads to the relation 
K = {2 — g) / {271 — 2 + g) . Supposing moreover that the parabolic solution extends almost over 
a whole period and that the correction brought about by the rounding off of the singularity 
around the maxima of the potential is negligible, we also have 

{Vm - F^) ~ - = '^-^ (26) 



Now, subtracting equation ( pS]) from equation (21), we find in the limit T — > 0: 

^(Fm - F„) = (2 - g){yM - Vm) + -TlogT (27) 

Combining equations ( P5| ) and (P?]), we find that for the previous equation to have a fixed 
point as T — > 0, (7 ^ 2 with negative corrections as: 



g)^^2'K5T\og^ (28) 

This result is independent of the way we calculate g, the correction to the surface tension. In 
particular, it shows that at zero temperature, the surface tension diverges as (L/a)^, where 
L is the size of the system. 

We can deduce an equation satisfied by JF', by plugging the ansatz for the derivatives of 
V in the fiow equation for V : 

dV' V -,2 V'" T V'" 

^ = {2-g- 2t,)V + 27r— ^ + W -^=^ + -— (29) 

i + v n + vy ^l + V 

— dV 

Close to the fixed point, we again suppose that to leading order in T, — — = in the above 

— / — // — 
equation. Plugging in the ansatz for V and V , and considering the leading term in T, we 

get to lowest order in T: 
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T + \ T + — T =0 (30) 



riu) r{u) 2T'{u 

We can show that necessarily a — b = 1 — a. Indeed, if 1 — a < a — 5, J-"" would be equal 
to zero while if 1 — a > a — 5, JF'(O) would be equal to zero, both alternatives being thus 
impossible. Hence, equation (^) can be rewritten as 



which yields after integration: 



^'t")log(|^) =2x«= (32) 



At this stage, we can note that the exponent relation 

2a- 5 = 1 (33) 

is independent of the scheme used to calculate of g (see Appendix). 

In the rest of this section, we calculate the exponents a and 7, and using g{T — 0) = 2, 
we also obtain JF'(O). These results now somewhat depend on the precise renormalization 
scheme we use to calculate the correction g to the surface tension. Replacing the derivatives 
of V by their expressions in terms of JF' and JF", in equation (K|), and changing variables from 

if to u = =a , we get to lowest order in T: 

From the exponent relation 2a — 6 = 1 derived previously, and the fact that g{T 0) is 
finite, we have another exponent relation a = 36, so that a = 3/5 and S = 1/5. We can also 
show that g{T 0) can be expressed in terms of JF'(O). From expression (^2]), one can see 
that JF' is a strictly increasing function on [0, 00] taking its values in [e, 00], and so we can 
change variables from JF' to its inverse function. Defining a new variable x as 

X = ^ , X (35) 
J^'(O) ^ ^ 



we have: 

u 

and 



^'^^^^^'\x\og{x)fl'' (36) 



27re 



du 1 log(x) / 27re \ 1/2 



dx 2 (xlog(x))V2 Vj^'(O) 
We can now express g{T 0) in terms of an integral over x from e to 00 as 
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Hence, using the fact that g{T — > 0) = 2, we finally find the constant JF'(O) ~ 1.38. 

5 Step energy as a function of temperature 

From physical considerations we know that below the roughnening temperature, the inter- 
face grows by forming terraces. An important quantity governing the kinetics of growth is 
therefore the step energy. The width of a step and its energy per unit length Ps can be 
obtained by comparing the elastic energy and the potential energy of a profile which 
changes by one period over the length ^. Requiring that these two energies are of the same 
order of magnitude leads to: 7/^^ ~ Vo where Vo is the amplitude of the periodic potential, or 

^ ~ a/—, and a step energy which scales as (3s oc Jv^. Since a step profile include Fourier 

V Vo 

modes such that < k < |A|, it is natural to use in the above equations the values of 7 
and Vo calculated for the length L = ae^ = ^. Since ^(L) ~ L/ ^yvo{L), one sees that this 
corresponds to stopping the renormalization procedure when Vo{L) ~ 1. We have integrated 
numerically the RG fiow, starting from 7 = 1 and from harmonic potentials of various am- 
plitudes Vo -C 1, and stopping for an arbitrary value V, chosen here to be Vc = 0.4. [] The 
resulting step energy as a function of temperature is plotted in Figure (H). For T close to Tr, 
one finds that ^ diverges as e^^^'^^~'^ , as it should since our RG fiow essentially boils down to 
the standard one [||]. For small temperatures, however, we find that f3s tends to a finite value 
with a linear slope in temperature. This slope is seen to decrease as the initial amplitude of 
the potential Vo increases. For Vo = 0.01, Ps decreases by ~ 30% when T increases from to 
0.25 Tr. This decrease falls to ~ 10% for Vo = 0.1. 

Experiments on Helium 4, on the other hand, have established that the step energy does 
only depend very weakly on temperature at small temperature, by not more than 5% when 
the temperatures varies from 0.05 Tr to 0.25 Tr. This suggests that the initial amplitude of 
the potential is of the same order as 70: in this case, the width of the step is of order a, and 
the bare parameters are not renormalized except possibly very close to Tr. The conclusion 
that experiments must be in the regime ~ 1 is in agreement with 0, where Vo is called 
tc (up to a numerical prefactor); fo/70 was estimated to be ~ 0.05. Since our RG fiow is 
different from the one obtained by Nozieres and Gallet, the values of the physical parameters 
obtained by a fit of our theory to the experiments will actually differ. 

6 Case of the contact line 

In this section, we repeat the previous analysis for the case of a contact line on a peri- 



odic substrate fl^, |14|. The roughness of a contact line on a disordered substrate, at zero 



^ Other values of Vc would not change the qualitative features reported below, provided Vc is not too large. 
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Figure 2: 13s/ \/v'o as a function of the rescaled temperature for three different values of the 
bare periodic potential Vo, with 7o = 1. 
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temperature, has been studied analytically and compared with the experimental predictions 
for the case of superfluid helium on a disordered cesium substrate, where the disorder arises 
from randomly distributed wettable heterogeneities which are oxydized areas of the substrate 
[ITSj p!6| . A physical realization of the theoretical situation we consider here could be achieved 
by preparing a substrate with equally spaced oxydized lines which would act as periodic 
pinning grooves. In this case the critical dimension is d = 1. We denote by $ the position of 
the line with respect to a mean position. The energy of the system is the sum of an elastic 
term and a potential term given by: 



= ^\k\mk)\^+ [\xvmx)) (39) 

where L is the length of the substrate and 7 the stiffness. 

The renormalization procedure is carried as before except that the propagator is now 

given by G(k) = — j — -. Moreover the renormalization of the stiffness now only comes from 
f3'y\k\ 

the scale change leading to the much simpler flow equation for 7: 

I- '«) 

— - — V - 2T 

Defining as before the rescaled parameters V and T with V = , ^, . , and T = — -r, the flow 

7A2|A| 7A2 

equations for V and T read: 

dV — V'^ T —ff 

and 

dT — 

During the flow, T flows to zero and the renormalized rescaled potential V develops shocks 
between which it has a parabolic shape. We characterize the singularities that develop around 
the maxima of V by the following scaling ansatz: 

l + F"(^) = TV^r(^^) (43) 

where ^'(0) > and ^(0) = 0. Putting u = g , this implies that for u ~ 1, 

V'{u) = -Te-ru + "f^^-^Hu) (44) 



and 



V{u) = Vm-T e T — + T e t g[u) (45) 
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with Q' = T . Plugging the previous expressions into the flow equation V : 

dV V" -,2 V'" T V" , , 

= -V + 2V =; + V =^ + =fj 46 

i + v {i + vy ^i + v 

dV' - 
and supposing that — — = 0, to leading order as T goes to zero, we have to leading order in 
_ di 

T: 

T e T + T e T + J- ^ T er = (47) 



r{u) r^{u) 2T'{u) 

Since JF'(O) > 0, we obtain a non trivial solution only if 



2a - (5 = 1 and A-B = B (48) 



and J-"' is again solution of equation (p2D. We can obtain the values of the parameters A 
and B by considering separately the singular part and the regular part of the renormalized 
rescaled potential V close to the flxed point. We expand V around one of its minima (p* as: 

Vi^) = V^ + ^icp-^r (49) 
and plug the resulting expression in the flow equation for V. This yields: 

^ =K. + -log(l + .) 

(50) 

dK 2k^ 



di 1 + K 

We note that the flxed point value for k is now flnite as T goes to zero and is given by k* = 1. 
Similarly the flow equation for V{0) = V m is: 

^=F„ + |log(T'e-#r(0)) (51) 

Combining equations (0) and (|5T|), the flow equation for the amplitude of V is given to 
leading order as T ^ by: 

^(Fa/ - Vm) = (Fm - F„) - I (52) 

Now, we expect that the singularity brings but a small correction to the parabolic part of the 
^* 

rescaled potential V , so that at the flxed point the amplitude of V is given by — , leading 

^ 

io A = 1. The value of the exponents a and 5 would require the analysis of subdominant 
terms. The conclusion of this section is that in the case of the contact line, the width of the 
singular region of the renormalized potential decreases exponentially with length scale: the 
potential quickly becomes a succession of matched parabolas. 
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7 Conclusion 

In this paper, we studied the problem of the thermal roughening transition using a frg 
formalism. We have shown that below the roughening temperature, the periodic potential 
on large length scales cannot be described by its lowest harmonic and that during the flow 
shocks are generated in the effective pinning potential. We expect that this result is more 
generally valid, and also holds in the case of a disordered pinning potential |jTl], [1^. By 
performing a resummation of our perturbation expansion, our results are in principle valid in 
the strong coupling regime, where the coupling constant is proportional to V j 7 (rather than 
V/T). Correspondingly, we stop the renormalisation procedure not when V{L) ~ T (as in 
ng), but rather when L reaches the size of the objects under investigation (for example the 
width of the steps). By comparing our numerical result with the experimental determination 
of the step energy of liquid Helium 4, we have concluded that the surface of Helium 4 crystals 
are such that the coupling to the lattice is of the same order of magnitude as the surface 
tension. This is in qualitative agreement with Balibar et al. 0, who estimate fo/70 ~ 0.05. 
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A Derivation of the flow equations 

In this appendix, we sketch the procedure to obtain the flow equations (|). We consider the 
partition function 

dm e-/5^('^(^)) (53) 

where Ti. is the hamiltonian given by (|l]). We spht the field $ into a fast moving and a 
slow moving component and average over the fast moving part. We can rewrite the partition 
function, up to a multiplicative constant, as: 



d[$<] e Uy^^^r (e J V A (54) 

where < ... >o represents the thermal average with respect to the gaussian weight: 



e J> (2^) (55) 
d'^k 

In the rest of this section we denote — r-j by dk. 

{2tiY 

A.l Renormalization of the periodic potential V 

We look for contributions to the potential V resulting from the above averaging, which are 
of the same form as the terms present in the hamiltonian before starting the renormalization 
procedure and which are of order di. These terms are represented by connected graphs and 
are all obtained by expanding the potential term with respect to $^ up to second order. 
There are only two ways of obtaining such graphs of order dl: 

• By contracting p two-legged terms — — d x — - — j V" \^ — - — j with 1 < p < oo. 

We must calculate: 



p\\ 2A2 



Jlld%V" 



p 



(56) 



which gives after averaging over the fast modes: 



Yldkjdk'j e'^'''^'''^^'''-'^''''^''>'^(^^>{ki)...^>{k'p 
i=i 



^\P / 1 \P f /^<(-r\\ r JL- ~ Jki{xi-X2)+...ikp{xp-xi) 
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In the above expression, the space dependence of V" is slowly varying, and since the integral 
is dominated by the region where the x^-'s are close to one another, we can with little error, 

treat these terms as approximately equal to V" — — J . After integrating over the rest of 
the ajj's and summing up over p, we are left with: 

- m I d'^x log(l + (58) 

d^x ^ — — y'y — X — / ^^^^ ^ two-legged 

terms I d'^x — V" — with < p < oo. To illustrate our method, we 
begin with p = 0. Expressing the fast modes in Fourier space, we have in discrete space: 



X y k k' 

which gives after averaging over the fast modes: 

X y A; ' ' 

In the above expression, the main contribution comes from the part x = y, while the part 

27r 

with X ^ y, has a phase which averages to almost zero. Using the fact that lAI = — , the 

a 

result is: 

Proceeding in a similar way for 1 < p < oo, we have 

f-f^Y ^ (p+l)(p + 2) /ly+Yay^+^)'^ 

^ V 2 ; (p + 2)! 2 \xy \l) 



j2 

Qikx+ik'y+i{ki+k{)xi+...+i{kp+k'p)xp (^$> (/^-^ $> (/^-'^ $> (^;^^) . . .$> (^^^) 

A/jfc'j/ci... Alp 

which yields after retaining the x = y part and averaging over the fast modes: 



(62) 
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J A — 1 A — 1 



(63) 



P _ _ gifc(x— xi)+iA;i(a;i— X2)...+ifcp(a;p— x) 
n ^^J^^ |^|2|L 12 12 



Treating the above expression as equation(|57D we are left with: 
Summing up over p, we finally obtain: 



2 y 7A2|A|2 + 7A2|A|2 ) 



(65) 



Taking into account the rescaling of the potential term, and supposing we are in the fiat 
phase so that $ and A are not rescaled, we obtain for d = 2: 

y,2 



7A2|A|2 

— V 

This fiow equation can be rewritten in terms of the rescaled parameters V = , ^, . ,^ and 

7A2|A|2 

T = ^ . Putting g = -^^we have: 
27r7A"^ 7 at 

A. 2 Renormalization of the surface tension 7 

The contributions to the gradient term are obtained from equations (0) and (§3|). 

• Consider first the contribution due to equation (^). Here, x is fixed since we have 
imposed that it should be equal to y. We can thus expand the V" terms with respect to this 
variable as: 

r (iiM) = V"{^) + ife - X) . Vt>(x)V""(q^) (68) 
The contribution of equation (|63D to the gradient term is given by: 
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n<m j=l 



P _ _ ^ikxi+iki(xi—X2)---+ikpXp 

n ^^^^^ |A;|2|A;i|2...|fcJ2 



(69) 



After integrating over Xj for i ^ m,n, we are left with: 

N 



2 ^ ^ 



n<m u=l 



^IM^ T/"p-2/'^lMV'"2/'^^) (70) 



dkdk'dk" 



2^nQy!^/|2^m-nQ^//|2^p+l- 



9 

Using the fact that ix^e^^^ = ^^e*'^^, and an integration by parts, we finally get after 

OK 

summing over n < m and over p: 



^/2/ $<(x) \-///2/ 'l'<(x) \ 
- ^di ^M^^yJ d'x ^ $<(x) ^ (V$<(x))2 (71) 



Since we are looking for the contribution to the gradient term, only the projection of the 
periodic function f on the zeroth harmonic counts. The contribution of equation 

(|63D to the elastic constant 7 is thus: 

,,,ij,,2.r /'..Q^ (72) 



d 



• Consider now equation (0). In order to obtain a term of the form (V$^(a;))^, we have 
to expand two V" terms. We proceed as follows: if we choose to expand y'($^(x„)/A) and 
l/"(<l><(xm)/A) with m < n, we perform the expansion with respect to (x^ + a;„)/2. The V" 
terms thus give: 
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Integrating over xj for j 7^ m, n, we get: 



2 P 7 V7A^/ 



m<n 



dkdk {Xjn Xn)jj 



|^|2)p-n+m(|^/|2). 



d d I 

Writing (x^ - x„)^, e'^''"^')^^™-"^") = _^_^e«(*:-*:')(a;m-a;„)^ g^^j performing the rest of the 

calculation as described above, we find that the contribution to the elastic term in terms of 
the rescaled parameters V and T is given by 



- / (75) 

finally leading, for = 2, to the renormalization of 7 given in the main text. If we had chosen 
another expansion scheme to obtain the contribution to the gradient term, for instance if we 
had expanded the terms with respect to the centre of mass, we would have obtained a 
somewhat different value for g. This would only affect the precise value of the exponents a 
and 5 obtained in the text, but not the qualitative features of the solution. 
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